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Starting from a microscopic definition of an alignment vector proportional to the polarization, 
we discuss the hydrodynamics of polar liquid crystals with local Coou-symmetry. The free energy 
for polar liquid crystals differs from that of nematic liquid crystals {Daoh) in that it contains terms 
violating the n — > — n symmetry. First we show that these .Z2-odd terms induce a general splay 
instability of a uniform polarized state in a range of parameters. Next we use the general Poisson- 
bracket formalism to derive the hydrodynamic equations of the system in the polarized state. The 
structure of the linear hydrodynamic modes confirms the existence of the splay instability. 
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I. INTRODUCTION 

In the past few decades, our understanding of complex 
materials has greatly benefited from the increasing so- 
phistication of technology and experimental designs. De- 
pending upon the various properties of the constituent 
molecules, such as geometric shape, electric and mag- 
netic moments, as well as chemical affinity, we are able 
to observe a variety of phases beyond the traditional crys- 
talline and liquid phases. Examples range from the tra- 
ditional nematic and smectic fiquid crystal phases to 
the recently discovered smectic blue phase or the ba- 
nana phases The possibility of ferroelectric order in 
liquid crystals has attracted much attention both from a 
fundamental and practical viewpoint. Among the com- 
monly observed liquid crystalline phases, only the chiral 
smectic-C* phase is known to be polar. Ferroelectricity 
in this phase was discovered in 1975 0] and the related 
phenomenon of antiferroelectricity in liquid crystals was 
discovered in 1989 Q. In these chiral liquid crystals, 
polarity is realized via the introduction of chiral carbon 
atoms into achiral molecules. Ferroelectric switching was 
also recently confirmed in the banana smectic phases Q 
consisting of achiral bend-core molecules. The packing 
structure in these banana phases induces a polar order 
in the layers of these liquid-crystal systems. 

The possibility of ferroelectric and antiferroelectric or- 
der in achiral polar liquid crystals composed of molecules 
carrying strong permanent electrical dipole moments 
has been a long-standing theoretical and experimen- 
tal question. Experimentally, ferroelectric order due 
to dipole-dipole interactions between rod-like aromatic 
copolyesters molecules has been reported jE S B S ^ 
by the Tokyo Tech group. In this case the nematic liq- 
uid crystal is biaxial and polarity was observed along 
both symmetry axes via second-harmonic generation. 
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mean-field) and numerical calcu- 
ITH IT^ Ir^ . as well as simula- 
have argued that there is no funda- 
mental reason forbidding the establishment of ferroelec- 
tric order in polar liquid crystals, assuming appropriate 
surface stabilization. 

The simplest ferroelectric liquid crystal of this type 
would be a uniaxial liquid of molecules carrying a per- 
manent electri c dip ole moment which lies parallel to the 
molecular axis [22| ■ In this case ferroelectricity would be 
due solely to dipolar interactions, with all dipoles point- 
ing in the same direction in the ordered state. In their 
classic textbook de Gennes and Prost argue that such 
dipolar molecules would always be asymmetric in shape, 
with head and tails of different size. This would naturally 
yield a splay deformation, making the liquid ferroelectric 
state unstable ^|. Theoretical work has indeed shown 
that any type of molecular asymmetry will suppress ferro- 
electric order in nematic liquid crystals, possibly replac- 
ing it with modulated polar phases [l7l |. In this paper 
we show that the instability of an ordered ferroelectric 
phase of uniaxial polar liquid crystals is more general 
than originally suggested by De Gennes and Prost. We 
demonstrate that in uniaxial polar liquid crystals a fer- 
roelectric phase is theoretically possible but only for a 
limited range of parameters. Outside this range the uni- 
formly polarized state is unstable to splay distortions, 
regardless of molecular shape. This may explain the dif- 
ficulties encountered in establishing experimentally the 
existence of a ferroelectric polar state. We employ a 
standard Poisson bracket formalism to derive the 

hydrodynamic equations for a liquid of uniaxial rod-like 
molecules carrying permanent dipole moments. From 
these equations we obtain the hydrodynamic modes in 
the uniformly polarized state and show that in a range 
of parameters the ordered state becomes unstable against 
the growth of long- wavelength splay fluctuations. This is 
consistent with an analysis of the free energy which in- 
dicates the presence of a mechanical instability for the 
same parameter range. Our derivation is based solely on 
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symmetry arguments and makes no assumption on the 
shape of the molecules. 

The familiar nematic liquid crystal phase is character- 
ized by a broken rotational symmetry due to macroscopic 
order of the rod-like molecules along a fixed direction in 
space. The broken symmetry direction is identified by a 
unit vector n, known as the director. The states of di- 
rector n and — n are indistinguishable. If the molecules 
carry a permanent electric dipole moment, in the ordered 
state there as many 'up' dipoles as 'down' dipoles and the 
system is not ferroelectric on large scales. The main dif- 
ference between a polar liquid crystal and the familiar 
nematic phase is the breaking of such a 2^2-symmetry 
(n — n). If the molecules carry a dipole moment, on 
average the dipoles are pointing in the same direction in 
the ordered state and the system can be ferroelectric on 
macroscopic scales. The order parameter characterizing 
the transition from the isotropic liquid phase to the po- 
larized phase is the alignment vector P, proportional to 
the mean polarization of the system. This should be con- 
trasted to the nematic order parameter Q^ , which is a 
rank-two tensor, also known as the alignment tensor. In 
uniaxial systems the mean values of both the alignment 
vector and the alignment tensor in an ordered phase can 
be expressed in terms of a unit vector pointing along the 
broken symmetry direction. 



P = Pon , 



(1) 
(2) 



The unit vectors n and n' are the directors in the polar- 
ized and nematic states, respectively. A macroscopic po- 
larized state is characterized by Pq ^ 0. From this it fol- 
lows that 5 7^ as well in the polarized state, and n' ~ n. 
The nematic phase has Pq = 0, but S' ^ 0. In the follow- 
ing we derive hydrodynamic equations for a polar liquid 
crystal. The hydrodynamic fields in the isotropic phase 
are simply the five conserved densities (number density, 
momentum density and energy density) . In the polarized 
phase we must add the two broken symmetry variables 
corresponding to the two independent components of the 
director n. In addition, we obtain the dynamical equa- 
tion for the full alignment vector P, which determines the 
order parameter for the polarized-isotropic phase transi- 
tion. Its magnitude is not a hydrodynamic variable, but 
it plays an important role in describing the dynamics near 
the transition, where its relaxation rate may become very 
slow. 

Our derivation starts in section II with the micro- 
scopic definitions of the relevant hydrodynamic fields. 
Using standard definitions from classical mechanics, in 
section III we obtain the coarse-grained Poisson-brackets 
for these macroscopic fields. Careful consideration is 
given to the derivation of the Poisson bracket between 
the alignment vector P and the momentum density g and 
that between the director n and g. A similar approach 
was taken in recent work on deriving the low-frequency 
hydrodynamics of nematic liquid crystals and nematic 



polymers [23, |27| . In section IV and V, we obtain the 
dynamical equation for the director n and the alignment 
vector P, as well as for the conserved densities. These 
equations contain new terms violating the 2^2-symmetry 
of the nematic phase. Finally, in section VI we exam- 
ine the hydrodynamic modes and demonstrate the splay 
instability induced by the new Z2-odd terms. 



II. DEFINITION OF MICROSCOPIC FIELD 
VARIABLES 

We start by defining the microscopic quantities of in- 
terest. The system is composed of N identical liquid 
crystal molecules, indexed by a. Each molecule consists 
of s atoms treated as point particles, with and the 
mass and the charge of the /i-th atom, respectively. We 
denote by r"^ the position of the fi-th atom in the a-th 
molecule. The position of the center of mass of the a-th 
molecule is given by 



(3) 



and y^, = 
that the molecules are neutral, 



is the molecular mass. We assume 



(4) 



but possess a no n- vanishing dipole moment d" , given by 



J2 9^1-"^ 



where is a unit vector along the molecular axis. 



and 



(5) 



(6) 



(7) 



The microscopic mass density, /5(r), momentum density, 
g(r), and dipole moment density, d(r) are defined in the 
usual way as 



P(r) = 

a.fi 

g(r) = ^p"'^,5(r-r"n 
d(r) = ^d"(5(r-r") 



(8) 
(9) 
(10) 



where p"'' is the momentum of the ^-th atom on the 
a-th molecule. The macroscopic mean variables describ- 
ing the dynamics of equilibrium fluctuations in the sys- 
tem are obtained from the microscopic ones after coarse- 
graining, p(r) = [p(r)]^, g(r) = [g(r)]^, d(r) = d(r) 
as described, for instance, in Refs. [3| and [2^. 
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The macroscopic coarse-grained fields, denoted by 
$(j(r,t), are hydrodynamic or quasi- hydrodynamical 
variables whose characteristic decay times are much 
longer than the underlying microscopic time scales and 
diverge in the limit of long wavelength. They are either 
conserved variables or variables associated with broken 
symmetries of the system. They evolve in time accord- 

ing to iEa 



dt 



(11) 



where Tab is the dissipative tensor, Ti. is the systems 
Hamiltonian, and y is the reactive term or nondissi- 
pative velocity given by 



K(r) = - J d\'Pab(v,v') 



sn 



(12) 



The Einstein summation convention on repeated indices 
is understood and 



Pab{r,r') = {$,(r),<i>,(r)} 
= -Pba(r',r) 



(13) 



denotes the Poisson bracket of the coarse-grained vari- 
ables. This is defined as 



with 



{$,(r),$,(r')} = {4,(r),|.b(r')} 



(14) 



do 



dA{r) dB{r') 
or - ^ op - ^ 



(15) 



The kinetic tensor F^i^ must be symmetric and 
9<i>^(r, t)/dt and STC/S^i, must have opposite signs under 
time reversal. The coarse-grained dynamical equations 
()11|1 can be rig orously derived from microscopic princi- 
ples liiliSlla- 

Polar order is described by the alignment vector P(r), 
defined as 



P(r) = ^"'^(^ 



(16) 



The alignment vector P(r) embodies the fiuctuations in 
the dipoles' orientation. We also write 



P(r) EE p(r)S'(r)n(r) . 



(17) 



The scalar field S{r) encodes information on the degree 
of local alignment of the dipoles, i.e., S{r) is the order 



parameter for the isotropic-ferroelectric nematic phase 
transition. The unit vector n(r) is the macroscopic 
director. As previously mentioned, the vectorial nature 
of P(r) is in sharp contrast to the tensorial character 
of the usual nematic order parameter, Qij(r). This 
mathematical distinction leads to new terms in the free 
energy of a polar fluid that are not allowed for a nematic 
liquid crystal. 



III. FREE ENERGY AND STABILITY OF THE 
POLARIZED PHASE 

Before proceeding with the derivation of the hydrody- 
namic equations for the polarized state, it is instructive 
to discuss the coarse-grained free energy of a polar fluid. 
This is given by 



F = Fp + Fn, 



Fp ^ f d\ sMll V • n(r) 
J Po 



F„. = 



l-Ksin X V X n)2 + . ..} , 



(18) 
(19) 

(20) 



where Sp{r) — p(r) — po is the deviation of the local den- 
sity from its equilibrium value, pQ. The contribution Fn 
is the free energy of a nematic liquid crystal which con- 
tains only terms symmetric with respect to n ^ — n, with 
Ki , K2 and the usual Frank constants for splay, twist 
and bend deformations, respectively. The contribution 
Fp contains additional terms that break that symmetry 
and are only allowed in a fluid that is macroscopically 
polar. In general a term ~ V • n(r) is also allowed in the 
free energy. This will give rise to a surface term that will 
favor a splay distortion of the molecules (corresponding 
to a non-zero V • n(r)). In the absence of surface sta- 
bilization, this would always destroy the polar order. In 
this work we assume that there is sufficient surface stabi- 
lization to suppress the surface term and to prevent the 
associated splay. 

In addition to a possible surface instability, it is ap- 
parent from the form of the free energy that a polar- 
ized liquid crystal can also exhibit a hulk instability to 
splay deformations. The latter is associated with the 
term ~ 5p{y)V ■ n(r) that cannot be eliminated by sur- 
face stabilization. The existence of this instability can 
be understood by rewriting the terms in the free energy 
involving density fluctuations and splay deformations of 
the director as 
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where V — {W±,dz). It is then clear that the effect of 
the couphng of density fluctuations to splay yield a down- 
ward renormalization of the splay elastic constant, Ki. 
For \B\ > \JK\C\ splay fluctuations becomes energet- 
ically favorable and will clearly destabilize the ordered 
state. The instability is suppressed if the fluid is incom- 
pressible since the renormalization of the Frank constant 
K\ vanishes when Ci — *■ oo. The stable ground state for 
\B\> \JK\C\ is expected to be characterized by director 
configurations that are spatially inhomogeneous in the 
X direction, perpendicular to the broken symmetry di- 
rection, with associated spatial structures in the density. 
The precise director configurations will depend on bound- 
ary effects, as well as on applied external fields (electric 
field in the case of ferroelectric order). The dynamics of 
the polarized state discussed below will of course reflect 
the existence of this instability, which will lead to the 
growth of the hydrodynamic mode associated with splay 
fluctuations for \B\> \/K\C\. 



In order to examine the possible relevance of the splay 
instability to experimental systems we need to estimate 
the various elastic constants. We consider a fluid of 
rods of linear size / interacting via a short-range repul- 
sive interaction of strength [/. The rods carry a per- 
manent electric dipole moment of magnitude d, which 
yields a dipolar coupling of strength Ud ~ kc(P'/r^ 
between two molecules at distance r, where is the 
Coloumb's law constant. The bulk compressional mod- 
ulus Ci has dimensions of energy density and its size is 
controlled by the repulsive part of the interaction, with 
Ci ~ U /l^ . The splay Frank coefficient has dimensions of 
energy /length and can be estimated as Ki ~ U/l. This 
yields an estimate for the critical value i?o — \/C\ K\ 
where the splay instability occurs: Bq ^ U/P. The elas- 
tic constant B that couples density and splay fluctuations 
has dimensions of energy/(length)^ and is controlled by 
interactions that favor alignment of the liquid crystal 
molecules (dipole-dipole interaction in the case of ferro- 
electric order). Clearly the elastic constant B would van- 
ish if the molecules are just rods as in the nematic case. 
We estimate B ~ Ud/I"^, where Ud is the strength of the 
dipolar interaction (for dipoles at a distance r ^ ^^'^ , 
we estimate Ud ^ kcCppo). The splay instability will 
then occur for Ud ~ U. Taking po ^ lO^^^m^^ and d ^ 
2 Debye we estimate Ud ^ 2.5 meV. This is smaller 
than U ~ O.leV, so we see that a stable ferroelectric 
phase of a polar liquid crystal is, in principle, possible. 
Of course, this assumes appropriate surface stabilization 
as discussed above. 



IV. POISSON BRACKETS FOR A POLAR 
LIQUID CRYSTAL 

We now return to the derivation of hydrodynamics 
via the Poisson-bracket method. Using Ea. H15() . it is 
straightforward to evaluate the Poisson-bracket relations 
between the various hydrodynamic fields. Our goal is to 
obtain equations describing the dynamics at long wave- 
lengths. We will therefore only keep terms of lower order 
in the gradients of the hydrodynamic fields and expand 
the 5-function as 

(5(r-r"^) = (5(r-r"- Ar"'^) 

= (5(r-r")- Ar^^Vfc(5(r~r") 
+0(V2) . 

(22) 

where Ar"^ = r"'^ — r". The required non- vanishing 
Poisson-brackets are: 

{p(r),5.(r')} - V.J(r-r')p(r'), (23) 
{5.(r),g,(r')} = g.(r')V,<5(r - r') 

-V:[5(r-r').9,(r')] , (24) 
{P.(r),g,(r')} = V, [P,(r)5(r - r')] 

-%Pfc(r)Vfe(5(r-r') 

+ \^Jk{Y)^k5{v~Y) , (25) 

where 

Ayfc = Y.^r^fDZ5{v-v^) (26) 

is a symmetric third-order tensor that depends on the de- 
gree of molecular alignment. We assume that its coarse- 
grained counterpart A^fe depends only on the alignment 
vector P(r). Quite generally, A^fe can then be written as 

PiP Pk 

Kjk = Ai ((5yPfc + SjkPi + SikPj) + A2 — , 

(27) 

where P{y) = n • P = pS. The two coefficients Ai and 
A2 are not independent. By contracting both sides of 
Eq. 1)2 7() with suitble combinations of 5ij and rik we ob- 
tain 

UkKk ^pS^ (4Ai + X2)pS , (28) 
n^njUkXijk = ^i'^a ■ nf S{r - r^) 

a 

= (3Ai+A2)p5. (29) 
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The two coefficients can then be expressed in terms of 
single microscopic quantity, A, given by 



A = 



p(r)5(r 



as 



Ai 
A2 



1~A 
5A- 3 



(30) 

(31) 
(32) 



Equation. (|27|l can then be rewritten in terms of A as 
1 - A 



A, 



ijk 



A njnfcFj(33) 



where Sfj is an operator that projects transverse to the 
director, 



(34) 



The director field itself has no microscopic definition in 
terms of the canonical coordinates and momenta of the 
individual atoms. Its Poisson-brackets must therefore be 
obtained indirectly from those of the alignment vector P. 
Using the definition of the director given in Eq. (|17|l and 
the chain rule of derivatives we can write 

{P,(r),5,(r')} = p(r)5(r){n.(r),5,(r')} 

+n.(r){p(r)5(r),5,(r')} • (35) 

Since ni{r) {^^(r), $p(r')} = 0, for any field <i>^, we ob- 
tain 



{n,;(r),.gj(r')} = 



p(r)5(r 



-C(r)m(r),5,(r')} (36) 



V. HYDRODYNAMIC EQUATIONS FOR THE 
POLARIZED PHASE 



In this section we derive the hydrodynamic equations 
in the polarized phase. The order parameter <S'(r) is as- 
sumed to be finite and constant in the following discus- 
sion. To calculate the reactive part of the hydrodynamic 
equations defined in Eq. H12I) in terms of the Poisson 
brackets, we need the coarse-grained Hamiltonian for the 
system 



n[p,g,n] = J d\ ^^ + F[pM 



(37) 



It consists of a kinetic part and a free energy F[p, n] that 
depends on the specific system in question and will be 
specified below. 

The nondissipative velocity for the mass-density equa- 
tion is immediately found to be 



VP 



-Vj5j(r) 



(38) 



The determination of the non-dissipative term for 
the momentum-density equation is straightforward albeit 
more tedious. Using Eqs. H23|) and H3()|l . we obtain 



-V 



.9»(i").9j(r) 



- p(r)V, 



SF 
6p{r) 



-Vfc 



.T SF 
SF 



Sni{v) 



Sji Xiki 



SF 



(39) 



p(r)5(r) Snj{r) 

Similarly for the director, we obtain using Ea. H36|l . 

U" = -VjV jUi + LOijTij + XSfi^Ukjnj , (40) 

where Vi = STi/Sgi is the velocity and Uij ~ i(ViWj + 
Vjfi) is the strain-rate tensor, while cuij — ^(Vifj — 
y jVi) is the vorticity. Only the part of the tensor Ay^ 
that is transverse to the director enters in the hydrody- 
namic equations. This is given by 



Sji{r)X.kiir) 



^p{r)S{r) 
X [Sljir)nk{r) - 



Sjkir)ndr)] 



(41) 



We now turn our attention to the dissipative part of 
Eq. 1)11(1 . Since the mass-conservation is exact, there 
is no dissipative term allowed in the dynamical equation 
for the mass-density. In the case of unbroken Galilean 
invariance, the time derivative of g can only couple to 
the gradients of velocity. As in the nematic case [23 , we 
introduce the tensor of viscosities 77^^; and the viscous 
stress tensor afj, with the following properties: 

'^jVijkl^kVl 



(42) 



where 



a5+ ae 



2 

-piSijUkk 



{riiUjk + njUik)nk 
+ PiiSijUkTiiUki +Ukkninj) . (43) 



To the lowest order in gradients, there is no difference 
in the expression of aYj between the polar and nematic 
cases. Likewise, as in a nematic liquid crystal, the time 
derivative dn/dt can only couple to a term of the form 



(44) 



where 7^^ is a rotational friction. 

Thus, collecting the reactive and dissipative terms, we 
arrive at the following hydrodynamic equations for a liq- 
uid crystal with polar oricntational order 



dp 
dt 



-V-g(r) 



(45) 



6 



dgj 
dt 



P 5p{v) 



5F 



1 + A. 



tSF 



1-A. 



5F 



Via. 



(46) 



(47) 



We note that the hydrodynamic equations for the polar- 
ized phase have the same form as those for a conventional 
nematic, as given for instance in [23|, although of course 
with a different microscopic expression for A. The only 
differences between the two sets of equations come from 
differences in the free energy. 

Using the form for the free energy given in Eqs. H2()|l 
and l|19|l . the hydrodynamic equations of a polar liquid 
crystal are given by the continuity equation 145|l and 



dtni{r) = ~Vj{r)\7jnt{r) +ujijnj + XSfi^UkjUj 
+7<5,^^V,p(r)-7*S§^, 



(48) 



VI. DYNAMICS OF THE ALIGNMENT 
VECTOR 

In the more general case where S'(r) is not a constant 
in Eq. p7l) , we need to consider the full order parameter 
of the alignment vector P. The hydrodynamic equations 
for the alignment vector P are 



d 



dt 



-V-g, (52) 

SF 

-P.(r)(V • v) - A,,,Vfci;,(r) - 7'— (53) 



dt 



5i(r) 



„ gi(r)g,(r) , _ SF 



6p{r) 



Pkir) 



SF 



SF 



SPj{v) 



V,cr, 



SP^ 



SF 



SP,{r) 



(54) 



The equation for the director can be obtained from 
Eq. (|53|l by using the relation 



dtni{v) 



1 



p(r)5(r 



■5l{v)d,Pu{v) . 



(55) 



p Po 



where 



SFr, 
Sp 



Sn 



B , , 1 + A 

--(V,7^,)(V,p)-^— 
Po ^Po 

1-A 



2po 



-p(r)V. 

1 + A 



SFr, 



2 

1-A 



Vfc 



rikS. 



(V,n 
tSF, 



BVk [SfjUkVjp] 



S_Fn 

Srij 



riiS. 



Via. 



(49) 



(50) 



5F 

-JL = -KiVjiV ■n) + K2{V x{V X n)) . 



-(7^2 - K-i)ns{V snk){V jUk) 
+ {K2-K3)Vb {nbn,{V.,n^)) 



(51) 



The terms in the hydrodynamic equations that are 
unique to a fluid with macroscopic polar order are those 
proportional to the elastic constant B. These are for- 
bidden by symmetry in the hydrodynamic equations of a 
nematic. 



VII. HYDRODYNAMIC MODES IN THE 
POLARIZED PHASE 

In the isotropic phase the only hydrodynamic fields are 
the four conserved densities: number density and three 
components of the momentum (energy fluctuations are 
not considered here). The corresponding four hydrody- 
namic modes are those of a conventional isotropic liquid: 
two propagating sound wave describing the decay of den- 
sity and longitudinal momentum fluctuations, and two 
diffusive shear modes, describing the decay of the two 
tranverse components of the momentum density. We now 
consider the long wavelength, low frequency dynamics in 
a uniformly polarized phase, characterized by uniform 
equilibrium values po of the density, Vq = of the flow 
velocity and no = z of the director field. In order to 
evaluate the hydrodynamic modes of the polarized state, 
we expand the hydrodynamic equations to linear order 
in the fluctuations of density, Sp = p — po, momentum 
Sg ~ Pov, and director 5n = n — z about their equilib- 
rium values (to lowest order z • = 0). We will ignore 
any variation in temperature in our consideration. We 
consider wave-like fluctuations whose space and time de- 
pendence is of the form cxp[i q • r — icut] and choose the 
coordinate system so that the wavevector q lies in the 
(xz) plane, i.e., q — {qx, qz) and qy = 0. We then project 
velocity and director fluctuations along the axes of an 
orthogonal coordinate system defined by the unit vec- 
tors (z,qa;,z X qx), with qx = qx/q- The equations for 
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velocity and director fluctuations transverse to both the 
direction of polarization z and Qx defined as 



(z X qx 



'q ' 



Sn"^ = (z X qx) • duq , 
decouple from the others and are given by 



1- A 



(56) 
(57) 



(58) 



.1 + A 



dtv^q = ~iyTm^^,y+^^^KTiq)q^q.6nl ,(59) 



where q = (fe/q, Qz/q) and 

Kriq) - [K2ql + K3ql]/q^ 



(60) 



"Tiq) = — + — J . (61 

The transverse eigenvalues are identical to those of ne- 
matic liquid crystals and are given by 



These modes are always diffusive and describe the cou- 
pled decay of vorticity and twist fluctuations. For qz = 
the transverse modes decouple and are simply given by 



lK2q^ , 

2/90 



(63) 
(64) 



where describes the diffusion of fluctuations in the 
local twist and controls the decay of fluctuations in 
the vorticity. 



= y ['^T(q) + 7^T(q) 



2 V Po 



(62) 



The remaining four hydrodynamic fluctuations, de- 
scribing the dynamics of density (5pq, velocity and di- 



rector components longitudinal to qx 

6n 



= qx • Vq and 

q — qx • Suq, and velocity component along the di- 
rection of polarization, are coupled. The linearized 
equations are given by 



dtSpq 



dtSfiq 



dtvl 



'iPmxVq - ipoqzVq , 

-jKLiqjq'^SUq + ij—qxSpd - i—^qzV^ 
Po ^ 



.1 + A 



-qxv„ , 



1 + A 



1 - A 



. 1 + A 

2po 



BqzqxSpq + i^-^^KL{q)qzq^5nl 



VL{ad<l «q - ^L^lxqzVq - iqx—Spq 

Po 

'^z{^)q'^Vq - v'^q^qzVq - iqz — Spq - ^T^BqlSpq - i^^^^KLiqjqxq'^Sn^ , 



1 - A 



(65) 
(66) 

(67) 

(68) 



where 



KLiq)^ [K^ql + K^ql]/q^ 



2(34 + as + ae qx 
Pi + P2 , 2^4 + a5 + ae 



Po 



4po 



(69) 



Q4 + Pl 9^ 2^4 + 05 + 0^6 ql .„„^ 

Po g"' 4po g"' 

ai + 04 + as + ae + Pi + 2p2 9^ 

f^4q) = 2 

Po q^ 



(71) 
(72) 



The corresponding hydrodynamic equations for a ne- 
matic liquid crystal are obtained from Eqs. H65I68|I by 
simply setting B = 0. For comparison we recall that 



there are four longitudinal hydrodynamic modes in a 
compressible nematic liquid crystal: two propagating 
sound waves and two diffusive modes describing the cou- 
pled decay of splay and velocity fluctuations 31]. The 
same four modes are present in a polarized liquid crys- 
tal, but their decay rates are modified in important ways 
by the new symmetry term breaking term proportional 
to B. 



We calculate the hydrodynamic frequencies up to o(q^) 
by solving the eigenvalue problem perturbatively. As ex- 
pected, we find two propagating and two purely diffusive 
modes. The two propagating eigenmodes are the sound 
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waves of the fluid and are given by 




;2 I ,, ~2 I r,,,/ -2-2 I -2 



Ci 



(73) 



To lowest order in the wavevector the propagation speed 



is proportional to 



As in isotropic liquids, it is de- 



termined entirely by the compressional modulus Ci of 
the fluid and receive corrections proportional to B. The 
coupling to polar order described by the elastic constant 
B increases the relaxation rate, but has no qualitative 
effect on the structure of the sound waves. In the incom- 
pressible limit, Ci tends to infinity and all contributions 
to the dynamics of the system from terms proportional 
to B vanish. This is evident in Eq. ()50|) as the -E2-odd 
term becomes a surface term. 

The eigenfrequencies of the two diffusive modes are 
given by 



Mff ^ 2 lK{c{) -I- /i(q) 1 
^ ^ 2 2^ 



PO 



(74) 



where 



M(q) = i^Mql + M^)ql ~ 

2q!4 -I- as -I- ae -I- -iaiqlql /4po • 



(75) 



(76) 



When K{q) = 0, the diffusive modes become unstable. 
This is precisely the splay instability discussed on Sec. 
Ill by examining the free energy of the system. 

In conclusion, we have derived the dynamical equations 
for the director n and the alignment vector P, start- 
ing from the microscopic expressions of the conserved 
quantities and carrying out rigorous coarse-graining via 
the Poisson-bracket approach. Our results are consis- 
tent with the equations for nematic liquid- crystals de- 
rived earlier using similar approaches |26l |2^ . Further- 
more, we have investigated the hydrodynamic modes of 
an ordered phase with macroscopic polarization. We have 
found that the polarity of the director field introduces a 
frustration in the system. As a result, a state of uniform 



local density and director orientation can become unsta- 
ble to the growth of long wavelength splay fluctuations 
in a region of parameters. The splay instability arises 
solely from the addition of 22-odd terms in the free en- 
ergy, which are generically allowed on the basis of sym- 
metry considerations in any polar liquid. We therefore 
believe that the existence of the instability is a general 
feature of any polar fluid, regardless of the origin of the 
polarity. 

It is our hope that the results reported in this paper 
will be useful in providing a rigorous justification of the 
various equations used previously in studying the ferro- 
electric phase of polar liquid crystals. It is our further 
hope that our results, presented here in the context of 
equilibrium physics, will also serve as a starting point 
in understanding the behavior of analogous systems in 
the realm of nonequilibrium physics where the polarity 
of these nematic phases is driven by active mechanisms. 
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Toner. W.K. and M.C.M. were supported by NSF Grants 
DMR-0219292 and DMR-0305407. 
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